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Abstract 

We consider a financial market with liquidity cost as in Qetin, Jarrow and Protter 
[3] where the supply function S e (s, v) depends on a parameter e > with S (s, v) = s 
corresponding to the perfect liquid situation. Using the PDE characterization of Qetin, 
Soner and Touzi [6] of the super-hedging cost of an option written on such a stock, we 
provide a Taylor expansion of the super-hedging cost in powers of e. In particular, we 
explicitly compute the first term in the expansion for a European Call option and give 
bounds for the order of the expansion for a European Digital Option. 
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1 Introduction 



The classical option pricing equation of Black & Scholes is derived under several simplifying 
assumptions. The "infinite" liquidity of the underlying stock process is one of them. In 
an attempt to understand the impact of liquidity, Qetin, Jarrow, Protter and collaborators 
[3, 4, 5] postulated the existence of a supply curve S(t, s, v) which is the price of a share of 
the stock when one wants to buy v shares at time t. In the Black & Scholes setting, this 
price function is taken to be independent of v corresponding to infinite amount of supply, 
hence infinite liquidity. In a recent paper, Qetin, Soner and Touzi [6] used this model and 
studied the liquidity premium in the price of an option written on such a stock with less 
than infinite liquidity. They characterized the option price by a nonlinear Black & Scholes 
equation, given in (2.3) below. In this pricing equation the liquidity manifests itself by 
means of a liquidity function £, which is given by 



£(t,s) :-- 



4^(M,0) 
av 



(t,s) € [0,T] x 



The liquidity function i measures the level of liquidity of the market. Namely, the larger £ 
is, the more liquid the market is. 

The main result of [6] is the characterization of the liquidity premium as the unique 
viscosity solution of a nonlinear Black-Scholes equation (2.3), which is very similar to 
the one derived by Barles and Soner [2]. This nonlinear equation can only be solved 
numerically as no explicit solutions are available. Motivated by this fact, in this paper we 
obtain rigorous asymptotic expansions for the liquidity premium. For vanilla options with 
sufficiently regular payoff, this expansion can be calculated explicitly giving further insight 
into the liquidity effects. 

As stated the chief objective of this paper is to analyze the large liquidity effect. Thus, we 
assume that the supply function depends on a small parameter e 

S £ {t, s, v) := S(t, s, ev), (t, s) G [0, T] xl + . 

Then, the corresponding liquidity function is given by 

£ £ (t, s) := -£(t, s), (t, s) G [0, T] x R+ . 

e 

Hence, as e tends to zero, the market becomes completely liquid. So we expect the price 
of an option V s to converge to the classical Black-Scholes price, v BS , and we are interested 
in expansions of the form 

ye = V BS + £V (1) + _ m + £ n v (n) + +0 ( £ «). 

Indeed, we prove this type of results and identify the functions i/ n ) in some cases. In 
particular, we show that 



Jt At< " Q ^ ^ SS ' 



" du. (1.1) 



This is exactly the liquidity premium of the standard Black-Scholes hedge 
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The paper is organized as follows. The problem is introduced in the next section and the 
approach is formally introduced in Section 3. Under a strong smoothness assumption, full 
expansion is obtained in Section 4. A quick convergence result is proved in Section 5. The 
Call option is studied in Section 6 and the Digital option in the final section. 

2 The general setting 

Let P) be a complete probability space endowed with a Brownian motion W with 

completed canonical filtration F = {Ft,t € [0, T]}, where T > is fixed maturity. The 
marginal price process St is defined by the stochastic differential equation 

^=a(t,S t )dW t , 

where a is assumed to be bounded, Lipschitz-continuous and uniformly elliptic. 

Given a continuous portfolio strategy Y with finite quadratic variation process (Y), the 

small time liquidation value of the portfolio is given by 

dzf = Y t dS t - [U £ (t, St)}' 1 d(Y) t = Y t dS t - e [U(t, St)]' 1 d(Y) t . 

The dependence of the process Z on its initial condition is suppressed for simplicity. 
Given a function g : R + — > R satisfying 

g is bounded from below and sup < oo, (2.1) 

s >o 1 + s 

the super-hedging cost is defined by 

V £ {t, s) := inf {z : Z £ t ' Y = z and Z £ / > g (S T ) F-as for some Y G A t , s ] , (2.2) 

where the time origin is removed to t and the initial condition for the price process is St = s. 
We refer to [6] for the precise definition of the set of admissible strategies At, s - 
This problem is similar to the super-replication problem studied extensively in [7, 8, 9, 
17, 18, 19, 20]. In the above setting, it is shown in Cetin, Soner and Touzi [6] that the 
value function of the super-hedging problem is the unique viscosity solution of the following 
nonlinear equation, 

-V t £ + H £ (t,s,V £ s ) = 0, on [0,T) x (0,oo), (2.3) 
satisfying the terminal condition V £ (T, .) = g and the growth condition 

-C < V £ {t,s) < C(l + s), (t, s) G [0, T] x R + , for some constant C > . (2.4) 
Here, H £ denotes the elliptic majorant of the first guess operator H £ : 

iT(M, 7 ) :=sup H £ (t,s,-f + P), 

/3>0 

H £ (t, s, 7 ) := -\s 2 a\t, S ) 7 - e[4£(t, a)]" W(t, s)^. 
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By direct calculation, it follows that 
H £ (t,s, 7 ) = -^s 2 a 2 (t,s) 



7+ 7 + 



2£(t, s) 



7+ 7 + 



For e = 0, both H £ , H e coincides with the following standard elliptic operator, 

H°(t, s, 7 ) = H°(t, s, 7 ) = -\s 2 a 2 {t, s) 7 , (t, s, 7) G [0, T]xl + x R. 

Hence, the equation (2.3) reduces to the linear Black-Scholes equation 

dv BS loo, x R<? 

-^-2*V(MK B / = o. 

We recall the well-known fact that its unique solution, v BS , is the Black-Scholes price, 



(2.5) 



BS 



(t,s) = E tjS [g(S T )\, (t,s)e[0,T]x 



where we used the notation E tjS = E[ • | St = s]. 



3 Formal calculations and Assumptions 

It is formally clear that as the market becomes more liquid, V s should converge to the 
Black-Scholes price v BS . Indeed, this is proved in Section 5. We are also interested in a 
Taylor expansion of V s in the parameter e, i.e., 

V £ (t, s) = v BS (t, s) + ev^\t, s) + e< 2 V(t, s) + ... + e n v( n \t, s) + o(e n ), (3.1) 

where o(e n ) is the standard notation, indicating that o(e n )/e n converges to zero as e tends 
to zero. 

Indeed, under sufficient regularity 



v (n) {t,s) 



1 d n V £ (t,s) 



de n 



e=0 



Thus, formally differentiate the equation (2.3) n-times with respect to e and then set e to 
zero. Using the above formal definition of y( n \ we arrive at, 





F n {t,s) 



—v. 



(n) 



Jt 2 
s 2 a 2 {t,s) 
U{t, s) 



-ls 2 a 2 (t,s)vif-F n (t,s), 



S) V 



(n-l-k) 



(t,s) 



(3.2) 
(3.3) 



k=0 



where we set := v BS . For all n > 1, the terminal data is v^ n \T, •) = 0, so that the 
Feymann-Kac formula yields 



,(»») 



n— 1 r „ 

(t, S) = J2 E t,s \ 



(u, S u )du 



(3.4) 



In particular, i/ 1 ) is given as in (1.1). 

The above calculations yield a rigorous proof when the pay-off is sufficiently regular. We 
will prove this in Section 4. On the other hand, for some discontinuous pay-offs the above 
functions may not be finite. For instance, for a digital option, t/ 1 ) = oo. Indeed, if we take 

g(s) := t s >K, a(t,s) = a and £(t,s)=£, 

we compute that 



( 1 i / a \ i o T-2u+t y 

I r T ( u _ £)g \<rs/T+u-2t Ui \K )^2 JT+u-2t) 

t (T-u)§(T + u-2t)§ 



( 1 ■ / s \ , o T-2u+t \ / 

fT e UvT+«-2t ln Un 2 VT +u -2tJ / ln(^) + cr T-2u + t 
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The first term above is actually +oo because of the non-integrability of (T — -u)~ 3 / 2 near 
T. 

In such cases, the expansion is not valid and a careful study of the behavior of V £ near the 
terminal data is needed. This will be done in Section 7. However, we first prove the full 
expansion in the "smooth" case. Then, in Section 6, we consider the Call option proving 
the expansion up to n = 2. Clearly, this later result extends to all Put options. Also, 
remarks on other payoffs and higher expansions are given in Remarks 6.2 and 6.1. 

4 Expansion for smooth pay-offs 

In this section, we prove the expansion under the assumption that there is a constant C so 
that 

-C < V W(t, s) < C(l + s), \(s 2 + l)v£\t, s)\<C, (4.1) 
\F n (t, s)\<C, V (t, s) G [0, T] x R+, n = 1, 2, . . . . 

Clearly, this is an implicit assumption on the pay-off g. Essentially, it holds for all smooth 
pay-offs growing at most linearly. In particular, (4.1) holds if a(t, s) = a, C(t, s) = I and if 
there exists a constant C so that 



-C <g{s) < C(l + s) 



an 

(» 2 + D^W 



<C, V s G R+, n = 2,3,. 



This is proved by using the homogenity of the Black-Scholes equation and differentiating 
it repeatedly. 

Following the techniques developed in the papers [13, 11, 14, 15, 16], for an integer n > 
we define, 

v ... (<i . )!= y(M)-saw.) i m 

where as before we set = v BS . 
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Theorem 4.1 Assume (4.1). Then, for every n = 1,2,..., there are constants C n and 
£o > so that for every e G (0, £q], and n = 1, 2, . . 



n-1 



.BS 



(t, s) < V £ (t, s) < v £ ' n (t, s) := E [e k v {k) (t, s)} + e n C n (T - t). (4.3) 

k=0 

In particular, as e 10, V s converges to the Black-Scholes price v BS uniformly on compact 
sets. Moreover, for every n > 1, V £ ' n converges to v^ n \ again uniformly on compact sets. 

Proof. Clearly, v BS < V s . We continue by proving the upper bound. Let v e ' n be as in 
(4.3) with a constant C n to be determined below. Using (3.2), we calculate that 

-v £ t ' n (t,s) + H £ (t,s,v^(t,s))>-v £ t ' n (t,s) + H £ (t,s,v £ s f(t,s)) 



= -vr-\sw s >r- 

n— 1 



2 2 

U(t, s) 



(v% n y 



fc=i 



es 2 a 2 
U(t, s) 



In view of (3.3), 



2 2 Ti— 1 

:r ' r ' fT " «f ) 2 - E ^ s )i = s ) + e " 



+i 



2^-2 



SO" 



4£(t, s) 



k=i 



U{t, s) 



9 £ (t,s) 



where g £ (t, s) is a quadratic function (t, s) for k < n and possibly powers of e. Hence 
by (4.1), there is a constant C n , 



n-1 



E &w,*)] 



k=i 



es 2 a 2 
U{t,s) 



< e n C n . 



Hence, we conclude that v £,n is a supersolution of (2.3). Moreover, by (4.1), — C < 
v £ ' n (t,s) < C(l + s). Then, by the comparison theorem for (2.3) (Theorem 6.1 of [6]), 
we conclude that V £ (t, s) < v £ ' n (t, s). 

In particular, this estimate implies the convergence of V s to v BS . To prove the convergence 
of V 6,n , we first observe that 



V s = E [e k v {n) {t,s)\ +e n V £ ' n . 



k=0 

Using the equations (2.3) and (3.2), we conclude that V £ ' n is a viscosity solution of 
-V £ ' n - \s 2 a 2 {t, s)V s T + F £ ' n (t, s, V S T) = 0, (t, s) G [0, T) xR + , 

where 



F £ ' n (t,s, 7 ) :-- 



H £ (t, s, v% n (t, s) + e n j) + \s 2 a 2 v £ s f + E ^ k F k (t, s) 



k=i 
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Tedious but a straightforward calculation shows that 

lim F £ ' n (t',s',j') = F n (t,s), 

where F n is as in (3.3). Then, by the classical stability results of viscosity solutions [1, 10, 
12], the Barles-Perthame semi-relaxed limits 

v {n \t,s) := liminf V £ ' n (t', a') and v^ n \t, s) := limsup V £ ' n {t' , a'), 
(t',s',s)->(t,s,0) (t',s',e)^(t,s,0) 

are, respectively, a viscosity supersolution and a subsolution of the equation (3.2) satisfied 
by v^ n \ Moreover it follows from (4.3) that 



(T, •) = y( n \T, •) = = v( n \T, •). 



We now use the comparison result for the linear partial differential equation (3.2), and 
conclude that > . Since 

v (n) (t, s) < lim inf V £ ' n (t, s) < lim sup V £ ' n (t, s) < (t, s) 

on [0, T] x R + , this proves that = = v^ n \ Hence, V £ ' n converges to the unique 
solution y( n \ uniformly on compact sets. 

□ 



5 A general convergence result 

In this section, we prove an easy convergence result under the following general assumption. 
We assume that 

cs 2 <£(t,s), (5.1) 

for some constant and 



Assumption 5.1 There is a decreasing sequence of smooth approximation g m > g of the 
pay-off g satisfying (4.1) with n = 1,2. Let Vm\ be the previously defined functions 
with pay-off g m . Then, F^(t, s) < c rn for some constant c m . 

This assumption is satisfied by all Lipschitz or for all bounded pay-offs. 

Theorem 5.1 Assume (2.1), (5.1) and that Assumption 5.1 holds true. Then, as the 
liquidity parameter goes to infinity, or equivalently as e \, 0, V s converges to the Black- 
Scholes price v BS . 

Proof. Let c m be as above and set 

u £ (t,s) :=v BS (t,s) + ec m (T-t). 
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As in the proof of Theorem 4.1, we can show that u £ is a super-solution of (2.3). Hence, 
V £ < u £ . Therefore, 

limsup V £ (t,s) < v BS {t,s). 



jBS, 



By (2.1), Vm S (t,s) converges to v BS (t,s). Since V s > v BS , this proves the convergence of 
V £ to v BS . 



□ 



6 First order expansion for convex payoffs 

One major limitation of our previous result is that the Call pay-off does not satisfy the 
Assumption (4.1). Therefore, in this section, we prove the first term in the Taylor expansion 
(3.1), i.e., 

V £ (t, s) = v BS (t, s) + evW(t, s) + o(e), (6.1) 

for convex payoffs satisfying weaker assumptions than (4.1). In particular, we will show 
that call options verify those assumptions. 

6.1 The general result 

In order to capitalize on the results we have already obtained for smooth payoffs, we will 
also consider a regularized version of our problem 



-V t e ' a + H £ (t, s, V S T) = 0, for (t, a) G [0, T)xl + , 
V £ > a (T,s)=g a (s), (6.2) 

where g a (s) = 4> a * g(s) with <j> a (-) := ^(^) and 4> is a positive, symmetric bump function 
on M, compactly supported in [—1, 1] and satisfying 



J 4>{u)di 



lu = 1. 

By convexity of g, for all a > we have g a > g, so that by monotony of our problem 



ye < ye,a^ 

Thus, since the main idea of our proof is to find a super-solution of (2.3), we see that it is 
enough to find a super-solution of (6.2). Let v BS ' a and v^' a , respectively, be the Black- 
Scholes price and the first-order expansion term for the regularized option. We now state 
our assumptions 

Assumption 6.1 (i) v BS + v BS ' a + + v < +oo. 
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(ii) As a tends to we have 

v BS > a (t,s) = v BS (t,s) + 0(a 2 ), 
v (1 ^ a (t,s) = v W {t,s) + o(l). 

(hi) There exists a constant c* independent of s, T — t and a and (v,f3) G [0, 1] x [1/2, 1] 
such that 1 < 2/3 + v < 2 and 

— (vW' a (t s)) 2 < - s \v BS ' a (t s)\ < - 

This assumption will be proved to be verified by Call options payoffs in subsection 6.2. 
Let V e > 1 be as (4.2), i.e. 

^'(M):= F ' (M '7 BS(M) . 

Theorem 6.1 Let Assumption 6.1 hold true and let a G (\, 2 f}+v ) • / or every (t, s) G 

[0, T] x R+ we have, 

V BS <V e < V BS ^ + + C*(T - t)/3+^ e 2-a(-+2/3) + C+(T _ t y € 3-2a(l+u)_ 

Moreover, V s — > v BS , V s ' 1 — > uniformly on compact sets, and (6.1) holds true. 

Proof. It is clear that V £ > v BS . To prove the reverse inequality, we start by following a 
technique similar to the one used in the proof of Theorem 4.1. Set 

v e,2 . = v BS,e" + ^(1),^ + _ t) f>+^ ,2-^+2/3) + C+(T _ t) u ^-2a{l+u) 

We calculate that for (t, s) G [0, T) x R+ 

- vf + £ e (t, s, vlf) > -vf + H £ (t, s, vf) 

ce 2 -*"^ s5 , e « (1)>e - 1 2 2 ej2 «V £;2 , 2 

= 7tt + -777 n v f ' — evi s a v' — — [v' 

(T-i) 1_/3_ ^ {T-t) l - v 1 1 2 ss u y ss ' 

„ e 2-a{v+2P) ,3-2a(l+i>) „2_2 „2_2 

i : : f ' n ,s n BS,e a (l),e a 2 



(T-i) 1 ^ 4^ v ss ' 2£ 

In view of Assumption 6.1 (hi), this quantity is always positive. We now analyze the terminal 
condition. In view of the conditions imposed on a, (3 and u 

v^ 2 (T,s) = v BS ^(T,s) = g £a (s). 

Hence, t> £ ' 2 is a super-solution of (6.2) and therefore of (2.3). Then, by the comparison 
theorem for (2.3) (proved in [6]), we conclude that V £ (t,s) < v £,2 (t,s). 

We now let e go to in the above inequalities. This proves that V s converges to v BS 
uniformly on compact sets. 

Finally, by Assumption 6.1(h) 

< V £,1 (t, S) < S) + O (>in{l-a(2£ +i ,),2-2a(l+,)}^ + {s 2 ^ 1 ) , 
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where it is clear with our conditions on a, j3 and v that the o(-) and O(-) above go to as 
e tends to 0. 

Using this estimate, we then prove the convergence of V s ' 1 exactly as in Theorem 4.1. □ 

Remark 6.1 Higher expansions can be proved similarly, provided that we extend Assump- 
tion 6.1 for n > 2. 

6.2 Expansion for the Call option 

In this section, we take 

g (s) = (s-K)+, a(t,s) = a, i(t,s) = l, 

and we verify that Assumptions 6.1(h) and 6.1(iii) are satisfied, since Assumption 6.1(i) is 
trivial. 

Straightforward but tedious calculations using the Feynman-Kac formula yield 



/ / 4>{u) exp ( — -d\(s, K + era, r) 2 ] du, 
2ttt J-i V 2 J 



where 



asy2 

8«vr J J_i J_i \Jv(2t - v) 



T = T — t, 



di(s, k, t) = — *-= \n(s/k) + \o\fi, 
ayt 2 

S(t, v, s, k) = — -j^=^= ln(s/fc) 



o\J2t — v 2 y/2r — v ' 

h a (T,v,s,k,x,y) =exp ^-(5(r, v, s, k) 2 + ^== (log (l + ^) + log (l + ^)) 

The following two propositions, whose proof is relagated to the appendix, ensure that 
Assumptions 6.1(h) and 6.1(iii) are satisfied 

Proposition 6.1 There exists a constant c* ; independent of s, r and a so that for all 
(v,P) G [0,1] x [1/2,1]: 

s \v BS ' a (t s)\ < - — (v W ' a (t s)) 2 < - 
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Proposition 6.2 As a tends to we have the following expansions 

yBS,^ s) = V BS { ^ s) + a 2___ J ^ mv 2 dv + O(o 4 ); 

v W< a (t, s) = vW(t, s) - a SRat ^ T J J <t>{x)<t>(y)\x ~ V\dxdy + o(q), 
where do(s,k,r) = — ^ln(s/A;) — \a^fr. 

Remark 6.2 It is not hard to show that the results of Propositions 6.1 and 6.2 hold for all 
convex linear combination of call or put options. However, we cannot use the above proof 
for, say, a call spread option whose payoff is neither convex nor concave. 

6.3 Numerical Experiments 

In order to have a better grasp of the liquidity effects, we also solved numerically (with 
simple finite difference methods) the PDE (2.3). We represent below the behaviour of the 
liquidity premium (that is to say V £ — v BS ) when the time to maturity t and the spot price 
vary 



Prime de liquidity V(t,s)-V R S(t,s) 




Figure 1: Call liquidity premium - T = 10, K = 15, a = 0.5, e = 0.1, I = 1 

In the above figure, the liquidity effect is strongly marked for ATM options and disapears 
quickly for ITM and OTM options. This was to be expected. Indeed, our calculations 
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showed that the liquidity effect is, for the first order, driven by the T of the call option 
(see (A.l)), which explodes for ATM options near maturity. Moreover, with our set of 
parameters, the first order correction is at most 0.06 for a BS price of 8.56, which means 
that the hedge against liquidity risk is not that expensive when the illiquidity is not too 
strong. 

We now compare the real liquidity premium with its first-order expansion term. 




Figure 2: Call first order liquidity premium - T = 10, K = 15, a = 0.5, e = 0.1, £ = 1 

A rapid examination of the above figure shows that the first order approximation remains 
excellent as long as we do not go too far from the maturity time T and we stay close to the 
money s = K. Otherwise, the first order overvalues the liquidity premium. 

7 Digital Option 

In this section, we analyze the specific example of a Digital option in the context of Black- 
Scholes model with constant liquidity parameter 

g(s):=l s >K, and a(t,s)=a, £(t,s)=£. 
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7.1 Theoretical bounds 



As pointed out earlier, for the Digital option, the first-order term that we obtained formally 
is equal to +00. Thus, the expansion (3.1) is no longer valid and our aim in this section is 
to find bounds for the first-order of the expansion. We start by approximating the option 
by a sequence of regularized call spreads. Then the original problem (2.3) is replaced by 

-V t e ' a + £ e (t, s, V S T) = 0, for (t, s) G [0, T) x R+, 
V e > a (T,s)=g a (s), (7.1) 

where g a (s) = 4> a * 9a(s) with g a (s) = ( s - K + 2a ) + ~( s - K + a ) + . 

Since <p a has compact support in [— a, a], notice that g a > g. Then, since the terminal 
condition is smooth, it follows from the comparison principle that 

V £ (t, s) < V e ' a (t, s), for (t, s, a) G [0, T}xR+x R* + . (7.2) 

With the same notations as in the previous section, we directly calculate using again the 
Feynman-Kac formula that 

1 f 1 

V g S ' a (t,s) = === / (j){u) f e -^(s,K+au-2a,T)2 _ ^d^K+au-a,^ ^ 

asa\/2iTT J -1 V / 

ft rl rl 



v m, a(t s) = _JL_ r f f ^)MM[^h^A dxdvdv 



where 



h a (T,v,s,K,x,y) = ^2 h a (r,v,s,K,x -i,y - j). 

1<M<2 

Then, we have the two following propositions which are proved exactly as in the call option 
case (since the functions involved here are essentially the same) 

Proposition 7.1 There exists a constant c* ; independent of s, r and a so that for all 
(v,P)e[0,i\ x [1/2,1] 

s \v BS > a (t s)\ < - ^-(v {l) ' a (t s)) 2 < - 

Proposition 7.2 As a tends to we have the following expansions: 

a -\d {s,K, T f 

v BS >"(t, s) = v BS (t, s) + 2 * e K(jV _ + 0{a% 

e -±d (s,K,rf rl rl 

V^>' a (t,s)=a l ^-^-j= J J 4,(x)4,(y)(\x-y-l\+\x-y+l\-2\x-y\)dxdy + o(a 1 ). 
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Define V e ^ c by 

r ....c (t , 8);= ^(M)-^(M) 



Theorem 7.1 Lei (/3,i/) G [1/2,1] x [0,1] be such that 7 := fpXl+l G (°> X ) anrf sei 
a := §(1 - 7). TTien for all (t, s) G [0, T] x R + , 



In particular, V s converges to v BS , uniformly on compact sets and 

0< liminf y £ ' 1 ' a (t',s / ,a) < limsup V £ ^' a {t' , s') < ^^+c*(T-t) ^) , 

(t',s',e)^(t,s,0) (*',«' ,e)->-(t,s,0) 2 K(T\/2lTT 

i.e. the order of the expansion is at least 2/5. 

Proof. It is clear that V s > v BS . To prove the reverse inequality, we start by following a 
technique similar to the one used in the proof of Theorem 6.1. Set 

v e,2 . = v BS,e" + ^(1),^ + ^ _ t) f) + *=± e 2-3a-a(„+2/3) + ^ _ t y e 3-2a(3+,) _ 

We proceed exactly as in Theorem 6.1 using Proposition 7.1. The result is 

-v £ t ' 2 (t, s) + i7 £ (t, s, v £ J(t, s)) > 0, for (t, s) G [0, T) x R + . 
We now analyze the terminal condition. Since 2/3 + ^ > 1, we have 

v £ ' 2 (T, s ) = u ss ' e °(r,s). 

Hence, f e ' 2 is a super-solution of (6.2) and therefore of (2.3). Then, by the comparison 

theorem for (2.3) (proved in [6]), we conclude that V £ (t, s) < v £ ' 2 (t, s). 

Then by Proposition 7.2 and the conditions imposed on a, j3 and v, we obtain easily the 

uniform convergence on compact sets of V £ to v BS by letting e go to 0. 

Now for the first order term, we would like to use our expansions and obtain a finite 

majorant for V £,1,c with the largest possible c. It is easy to argue that c = a is the best 

choice possible. This, in turn, imposes the following condition 



a < min 



1 



2'4 + 2/3 + i/7 + 2i/J 4 + 2/3 + z/ 



Now it follows that, for all 7 > small enough, there are /3 and v satisfying our conditions 
so that 4+2 2 3+i/ = |(l — 7). It suffices then to take the liminf and limsup in the inequality 
to prove the result. □ 



14 



7.2 Numerical results 



The digital option liquidity premium In this section, we provide numerical results 
for the case of the Digital option. As in the section 6.3 the PDE (2.3) is solved with finite 
difference method. We represent below the behaviour of the liquidity premium when the 
time to maturity t and the spot price vary 

Prime de liquidite V BS+|(a|pha) (t,s)-V BS (t,s)) 

05 *,-•••"""] 




Figure 3: Digital liquidity premium - T = 10, K = 25, a = 0.5, e = 0.1, I = 1 

Qualitatively, the liquidity premium behaves as in the Call case. However, as expected the 
effects of illiquidity are even stronger for ATM options near maturity, since the T of a digital 
option explodes faster. Moreover, with our set of parameters, the first order correction to 
the price is at most 0.04 for a BS price of 0.21, which means that the hedge against liquidity 
risk is much more expensive in the case of a digital option, for a same level of liquidity in 
the market. 

Numerical confirmation of the expansion order We represent below the liquidity 
premium for a fixed value of the spot when the parameter e varies with a logarithmic scale. 
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Figure 4: log (V £ - v BS ) - T = 1, K = 25, a = 15, <r = 0.5, e = 0.1, * = 1 

For small values of e we observe the expected linear behaviour of log (V s — v BS ) . The slope 
of the above curve is roughly equal to 1/2 (the exact value here is 0.54), which is close to 
our minimal value of 2/5. The numerical results suggest that the true expansion order lies 
in the interval [2/5, 1/2]. 

It is also important to realize the financial implications of our results. We just have high- 
lighted the fact that the first order effect exhibits a phase transition for discontinuous 
payoff, in the sense that derivative securities of the type of digital options induce a cost of 
illiquidity which vanishes at a significantly slower rate than the continuous payoff case. This 
means that derivative with discontinuous payoff are more rapidly affected by the illiquidity 
cost. 

Acknowledgements The authors whish to thank Reda Chhaibi for letting them use his 
Matlab code for the numerical resolution of the PDE (2.3). 
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A Technical Proofs 



Proof. [Proof of Proposition 6.1] We start by proving the inequality for vf s S,a . By dom- 
inated convergence, it is clear that svf s ,a goes to when s approaches or +oo. Hence 
for a 0, it also converges to when r tends to 0. Thus svf s S,a is less than a constant C a 
independent of s and r. However, when a tends to zero, we obtain the classical expression 
of the T of a call option 

-|dl( S ,^,r) 2 

vf s S (t,s) = -=^, (A.l) 

sa\/2irT 

which is known to explode only when s = K and r — > 0. Therefore, to understand the 
dependence in a of C Q , we only have to study the behaviour of svf s S,a when s = K and 
when both a and r go to 0. 

Let us therefore take a = e a and r = e b with a and 6 strictly positive numbers. For all 
(3 G [1/2, 1] we have 

,(6/2-a)(l-20) ,1 i(ae b ' 2 g^ 2 i J 1+ e a u\\ 2 

crV27r 7-1 



Therefore, if a < 6/2 (i.e. if r goes to faster than a) the quantity above always goes to 
when e — > due to the exponential term. If a > 6/2, the exponential term goes to 1, but 
since (3 £ [1/2, 1] the above expression has always a finite limit. Hence the inequality for 



BS,a 
SV SS . 



A change of variable and direct calculations imply that, for all v £ [0, 1], we have 

T^ a ^SV$-°(t, S) = /' / (A . 2) 

8«vrs 7 7(-i,i)2 vM2 - 



where 



^(r, v, s, K, x,y) _ n , f ^ ^ _ log (l + f) + log (l + - 



2 



h a (T,v,s,K,x,y) ~ y ' oV T ( 2 - u) / 

V CTVr(2-t;) / 

Using the same arguments as in the proof of the previous inequality, we can show again 
that the only problem corresponds to the case where s = K and a and r go to 0. Using 
the same notations, we have 
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x exp ^— log 1 + — log 1 + 



x exp 



a 2 (2-v) 6 V K J \ K 
e~ b (l og (l + ^)-log(l + f)) 2 \ 
2a 2 v(2 - v) ) 



h £ a(e\v,s,s,x,y) _ ( oei{l-2v) _ ft log (l + ^) + log (l + ^ 



2 



h £ a(e b ,v,s,s,x,y) \ y/2 - v a\/2 - v 

'ae b Hl-2v) b log(l + ^)+log(l + ^ y , . 

— ; ; + g~ fe v ^ cta/2 - vs?. 

V2 - v aV2 - v 

Therefore, if a < 6/2, /i £ a always goes to 0. Otherwise, the integral has a finite limite but 
since v £ [0, 1] and a > 6/2, the expression in (A. 2) has a finite limit. This proves the 
second inequality. □ 



Proof. [Proof of Proposition 6.2] The first result is straightforward and only uses the 
fact that the function <p is symmetric, which allows us to get rid off the odd terms in the 
expansion. For the second one, we directly calculate that 



f r rl rl } 

Jo 7-1 7-1 8^V(2r - v) 



v 



c 2 a 2 (x-y) 2 „(„2\ 

fr rl rl ~ b 4^^,(1- £) +°( a > § 

+ a Jo LL mm 8,iKaVv(2r-v) + 



t r l rl 



+ a L LL mm 16vrliW>(2r-,)3/ 2 

( r rl rl fP- J^^ ) \ 

+ o\a / / <j>{x)<j>{y) = Gtedycfo 

1 jo J-iJ-i 8tti\/v(2t — v) 



where we suppressed the arguments of the functions v^' a and 5 for notational simplicity. 
Note that all the above integrals are well-defined and finite. Then using dominated conver- 
gence and the fact that (f> is symmetric, it is easy to show that 
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V 



W> a = / / (f>(x)(/)(y) , dxdydv 

o 7-1 7-1 v / 87rfo(2r - w ) 

T fl fl „-5 2 i 



+ a l L L 8*£KaM2r -V) {X + v)dxdydV 
+ a L LL^ X)<KV) 16^V^(2r - W 2 dxdydV + °( a ^ 



R 2 a 2 (x-y) 2 , ( 2\ 

-T fl fl ~ S -«V,(l-») +(, ( a ) 



/■T fl fl e !K*<j*v(l-^) - ' 

= / / / 4>(x)(l>(y) - / =^ dxdydv + o(a). 

Jo J-i J-i 8itI\Jv(2t — v) 

Now the first term in the expansion above goes clearly to as a tends to 0. Then we 
have 



- «W = f Z 1 Z 1 ' e-^^» +0(a2) - 1 dxdydv + (a). 

7o 7-1 7-i 8it£Jv(2t - v) \ 



Using the change of variable u = 2Ku% ' ^ e ^ rs * term above can be rewritten as 



. . ~ S(-r-„-\i +°( a2 ) 



r +QO f 1 f 1 e - 5(T "4K 2 Jl 2 < S > K)2 (f){x)cp(y)\x -y\e ^Zrrtw -\ 



8-KlKo J a\x-y\ y_ 1 y_ 1 l~ a 2 (x-y) 2 V? 

■2K<?^F \ AT - AK 2„2,„2 



A simple application of the dominated convergence and Fubini theorems shows that the 
above integral (without the a factor) has a finite limit as a approaches and is given by 

e -\d a (s,K,T) 2 rl rl r+oo e -u 2 _ ^ 

Since the last integral is equal to y/n, we obtain the second expansion. □ 
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